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ABSTRACT 

We construct merger trees for galaxies identified in a cosmological hydrodynamical 
simulation and use them to characterize predicted merger rates as a function of red- 
shift, galaxy mass, and merger mass ratio. At z — 0.3, we find a mean rate of 0.054 
mergers per galaxy per Gyr above a 1:2 mass ratio threshold for massive galaxies 
(baryonic mass above 6.4 x 1O 1O M0), but only 0.018 Gyr" 1 for lower mass galaxies. 
The mass ratio distribution is oc i?mcrg f° r the massive galaxy sample, so high mass 
mergers dominate the total merger growth rate. The predicted rates increase rapidly 
with increasing redshift, and they agree reasonably well with observational estimates. 
A substantial fraction of galaxies do not experience any resolved mergers during the 
course of the simulation, and even for the high mass sample only 50% of galaxies 
experience a greater than 1 : 4 merger since z = 1. Typical galaxies thus have fairly 
quiescent merger histories. 

We assign bulge-to-disk ratios to simulated galaxies by assuming that mergers 
above a mass ratio threshold i? ma jor convert stellar disks into spheroids. With i? ma jor 
of 1 : 4, we obtain a fairly good match to the observed dependence of early-type 
fraction on galaxy mass. However, the predicted fraction of truly bulge-dominated 
systems (/bulge > 0.8) is small, and producing a substantial population of bulge- 
dominated galaxies may require a mechanism that shuts off gas accretion at late times 
and/or additional processes (besides major mergers) for producing bulges. 
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1 INTRODUCTION 

It has long been suggested that the spheroidal compo- 
nents of galaxies are a result of galaxy mergers (e.g., 
iToomre fc Toomrell973) . and that a galaxy's merger history 
therefore plays a central role in determining its Hubble type. 
Numerical simulations of disk galaxy mergers produce sys- 



ies (e.g.. Barnes 1988; 


Hernauisdll992: Barnes & Hernauist 


1996 


Naab et alj Il999l: ISDrineeJ 2000: 


Bendo & Barnes 


200C 


Cretton et alJl200lllNaab & Burkert 2003: Cox et al 


2004 


Springel & Hcrnauist 2004), includin 


g the telltale low 



surface brightness shells and tidal features that are of- 
ten revealed by deep imaging (e. g.. iMalin fc CarterlHiiS 
ISchweizer et al]|l99Ctlvan" Dokkum 2005). However, present 
simulations do not have the dynamic range needed to track 
the morphological evolution of individual galaxies while 
simultaneously modeling representative cosmological vol- 
umes, so it is unclear whether the merger hypothesis can 
explain the observed frequency and environment dependence 



of Hubble types. Semi-analytic models of galaxy formation 
have provided the main avenue for progress on this ques- 
tion, using simplified recipes to c onvert disk material in 
stellar bulges following mergers fe.g.lKauffmann et alj 1993; 
Kauffmann 1996; Baugh i e^alJ^^96jJSomer^ille_& Primack 
19991: ICole et alJl2000t iBell et al.ll2003al: lHatton et alJl200a 
N^gashitn^^^^shil2004ll . These models have achieved rea- 
sonable success in reproducing observed galaxy distribu- 
tions, but they rest on an approximate treatment of dy- 
namical friction to compute galaxy merger rates from dark 
halo merger rates, in addition to the assumptions of the 
merger recipes themselves. One recent semi-analytic study 
that tracks "galaxies" as sub-halos in high resolution N-body 
simulations yi elds results significantly different from those of 
other models (iKang et al]|2004lh In addition, the effects of 
dynamical friction and tidal disruption on sub-halos will be 
altered by the presence of tightly bound gas and stars within 
them. 

In this paper, we compute the merger histories of galax- 



2 Mailer et al. 



ies formed in a smoothed particle hydrodynamics (SPH) cos- 
mological simulation. While this simulation does not have 
the resolution required for reliable morphological analysis, 
it incorporates all of the physics that should have a major 
effect on merger rates, and we can apply semi-analytic style 
recipes to compute resulting distributions of Hubble types. 
We can also address, at least in part, the related issues of 
survival and angular momentum of galaxy disks. 

Stellar disks are fragile, and one important question 
for cosmological models is whether the y predict a low 
enou gh merger rates to preserve thin disks llToth fc Ostrikeil 
1992). In addition, disk galaxies forming in cosmologi- 
cal SPH simulations frequently lose angular momentum 
in late-time merger events, and the simulated galaxies 
have, on average, less angular m omentum than observed 
disks iNavarro fc SteinmetdHoOttl . Recent simulations have 
produced some examples of disk galaxies with realis- 
tic angular momentum properties iGovernato et all 120041: 
iRobertson et alJl2004h . but these results are strongly de- 
pendant on the galaxy's merger history. Therefore, with- 
out knowing the distribution of galaxy merger histories 
we cannot ascertain how often such realistic disks are cre- 
ated. It has also become clear that a halo's merging his- 
tory is related to its a ngular momentum iMaller ct al. 2002: 
IVitvitska et al]l2002ft . and the way that baryons trace that 
merging h istory may play an im portant role in determining 
disk sizes jMaller fc Dekelll2002l) . 

Our results for the merger statistics of SPH galaxies 
complement analy tic and numerical calculations of dark halo 
merg e r stati stics iKajjJfjnanii_fc_Wliijjll99j; |Lacev_fc_Coli 
1993 . Il994t ISomeryille fc Kolatti Il999t ISomerville et al. 
2000; Wechsler et al. || 2002r). They signif icantly improve on 
the earlier study of lt^rrali et all i2002T) by using a higher 
resolution simulation and, equally important, by tracking 
merger histories of individual systems rather than globally 
averaged merger and accretion rates. We describe the sim- 
ulation and our method of constructing merger trees in 
We present merger statistics in and in |J31 we use the 
merger histories to compute Hubble type distributions. We 
summarize our results and discuss directions for future work 
in 



BUILDING MERGER TREES IN THE 
SIMULATION 



The simulation was evolved using PTreeSPH (Davee^d 



1997), a parallel version of the SPH code of [Katz et al 
(1996), for a flat fl m = 0.4 cosmology with erg — 0.8, a 
Hubble constant Ho = 100/ikm s _1 Mpc _1 with h = 0.65, a 
baryon content Qi, = 0.047, and a spectral index n = 0.93. 
Our choices of as,Ho,n, and Qi, are reasonably close to the 
latest cosmological par ameters estimated from the CMB and 
large scale structure iSpergel et alJ 120031 : lEisenstein et alJ 
2005), while our value of Cl m is higher by about 1.5<r. The 
box is 22.22/i -1 Mpc on a side, with 128 3 dark matter parti- 
cles and 128 3 gas particles. The limited size of the simulation 
could have a significant impact on our statistics — in partic- 
ular there is only one cluster mass halo in our simulation 
volume, and that halo is itself anomalously large for a box 
of this size. The dark matter particle mass is 7.88 x 10 s M@, 
and the gravitational softening is a 5/i -1 comoving kpc cu- 
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Figure 1. The upper panel shows the cumulative galaxy mass 
function at z = for o ur simulat i on and , for comparison, the 
Schechter function fit of iBell et aT] i2003bl) . The y-axis is in co- 
moving units with h = 0.7 as we use throughout this paper. The 
simulations produce far too many low and high mass galaxies, 
while galaxies around the bend in the Schechter function are a 
factor of 2-3 too massive. The bottom panel shows the factor by 
which the simulated galaxies masses should be divided in order to 
agree with the observations. For the galaxy masses in which we 
are interested, 2 X 1O 1O M < A/ gal < 6 X 10 n M Q , the correction 
is roughly a factor of 2.75 as shown by the dotted line, and it does 
not depend strongly on galaxy mass. 



bic spline, roughly equivalent to a Plummer force softening 
of 3.5/i _1 comoving kpc. The gas particles start off with 
masses of 1.06 x 1O 8 M0, though this changes as star parti- 
cles are created. This simulation has p reviously been used 
for studying quasar absorption systems_j^^^^M I^S), the 
X-ray properties of galaxy groups JPave^tHdT 20021b the 
correlation bet ween Lyman Break G alaxies and Lya for- 
est absorption (Koll meier et alJ |2Q05j), and the histories and 
mech anisms of gas accretion (iMurali et all2002llKeres et alJ 
2005). Because the uncertainty in the Hubble parameter has 
been significantly reduced with recent observations and is 
now known to be quite close to h — 0.7, we scale all observ- 
able quantities to the values appropriate for h = 0.7 and do 
not quote explicit h dependences. (For example, we take the 
box volume to be 22.22/0.7 = 31.75 Mpc on a side.) 

We analyze 227 output time steps between z — 8 
and z = 0, spaced to give a time resolution better than 
150 Myr. At each step, galaxies are identified using the 



Spline Kernel Interpolative DENMAX (SKID ) 1 algorithm 



(S 

Gelb fc Bertschineerl Il99l iKatz et all Il996l) . This algo- 
rithm consists of five basic steps: (1) the smoothed baryonic 
density field is determined; (2) baryonic particles are moved 
towards higher density along the initial gradient of the bary- 
onic density field; (3) the initial group is defined to be the 
set of particles that aggregate at a particular density peak; 

(4) initial groups that are close together are then linked; 

(5) particles that do not satisfy a negative energy binding 
criterion relative to the group's center of mass are removed 



1 We use the implementation of SKID by J. 
Stadel and T. Quin, which is publicly available at 
http: / /www-hpcc. astro. washington.edu/tools/skid. html 
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Sample 


Number 




Mass Range 


R Mag. Range 


Min. i?merg 


Total Number 
Main Branch 


of Mergers 
All Branches 1:8 


2 res 

1:4 


1:2 


High 


166 


6.4 


- 60 X 10 10 M Q 


-18.3 < M R < -19.4 


0.125 (1:8) 


466 


595 0.5 


1.0 


2.6 


Medium 


167 


3.3 


- 6.1 x 10 10 M Q 


-17.6 < M R < -18.3 


0.250 (1:4) 


132 


137 


0.6 


1.0, 


Low 


165 


2.0 


- 3.3 x 1O 1O M 


-16.8 < Mr < -17.6 


0.500 (1:2) 


65 


67 




0.5 


All 


507 


2.0 


- 290 x 1O 1O M 


M R < -16.8 




752 


988 







Table 1. The properties of galaxies in our three mass bins are shown as well as the properties of the sample as a whole. "All" includes 
9 galaxies more massive then the upper mass cutoff of the high mass sample t hat are only used in § 13.21 The i?-band magnitude, 
corresponding to the same number density of objects in the B lan tern et alj i2003t) luminosity function is given for each mass bin. The 
minimum parent mass ratio for which the sample is complete and the total number of resolved mergers both along the main branch 
and for all branches are also given. For the medium and low mass samples essentially all resolved mergers are along the main branch. 
In the high mass sample about 30% of mergers occur off the main branch. We will focus on main branch mergers here as they are most 
important in determining the Hubble type of the z = galaxy. 



from the group. These steps are applied to all star parti- 
cles and to gas particles with temperatures T < 3 x 10 4 K 
and densities p g as/p g as > 1000. The resulting groups of cold 
gas and stars have masses an d sizes similar to the luminou s 
regions of observed galaxies jKatzlll992t iKatz et al.l fl996). 
The smallest systems reliably identified, as determined by 
comparison to higher resolution simulations of smaller vol- 
umes, have 64 particles, or a mass of 6.8 x 1O 9 M . However, 
for our merger study we restrict ourselves to galaxies with 
more than three times this mass, or A/ ga i > 2 x 1O 1O M0 
at z = 0, since we also need to identify their smaller mass 
progenitors. 

Previous merger trees that we are aware of have been 
built by identifying the progenitors (sometimes only the 
most massive progenitor) of each dark matter halo at time 
t\ in a previous time step fa, then repeating this proce- 
dure for all time steps |Lacev fc Colelll994l IWechsler et all 
| 2002l: IZhao et al.||2003|: iLin et alJl2003t iPeirani et alJl2004f 
ISheth fc Tormenl l2004t iTormen et all 12004) . However, we 
have found a number of complications when trying to apply 
this algorithm to simulated galaxies. The main difficulty is in 
reliably defining a galaxy in each time step. As is sometimes 
the case with optical observations, it can be difficult to de- 
termine if merging systems should be identified as one or two 
galaxies, and on rare occasions a galaxy is not identified at 
all in one time step. If a galaxy's progenitor is missed in one 
time step, then the above algorithm will truncate the tree 
and find that a galaxy suddenly appears in the simulation. 
More commonly, we find that two galaxies have the same 
progenitor, so that galaxies appear to split and then re-form 
into one galaxy (see the top panel of Figure|2J. A third com- 
plication is that particles can be stripped from both merger 
progenitors during the merging process but later join the 
merger remnant. Whether these particles should be counted 
as accreted mass or as part of the merger is again a matter 
of definition. These complications affect only a small minor- 
ity of galaxies, and the application of SKID to the star and 
cold gas distributions generally yields robust galaxy identi- 
fications. However, with 227 time steps, even an algorithm 
that is reliable 99.9% of the time will encounter problems 
for one out of five merger trees. 

We have developed a modified merger tree algorithm to 
address these issues. We distinguish between a SKID group, 
which is the output from running SKID on the simulation, 
and a galaxy, which we define in terms of our merger tree. 



Usually the two are identical, but occasionally a galaxy can 
be made up of two or more SKID groups, or contain parti- 
cles not in any SKID group, or both. We identify all SKID 
groups at z = 0, trace each group's members back through 
all previous time steps, and identify all progenitor SKID 
groups containing at least 16 of these members. 2 If a par- 
ticle is not part of a SKID group in a particular time step 
but is part of a SKID group at an earlier time step, then we 
assign it to a "virtual" SKID group. These virtual assign- 
ments allow a SKID group to have multiple descendants, at 
least temporarily, usually during and after mergers. 

To build the merger tree of a SKID group identified 
at 2 = (or another redshift chosen for the analysis), we 
trace its progenitors back in time but always combine all 
of the descendants (real and virtual) of a progenitor SKID 
group into a single object. This approach deals with the 
problem of temporarily split off particles or divided galax- 
ies. We do not count particles that are not in "real" SKID 
groups at 2 = 0. Occasionally the two descendants of a sin- 
gle galaxy never rejoin at a later time. When this occurs 
near z = 0, we interpret it as an interrupted merger process 
and anticipate that the two descendants would rejoin if the 
simulation were run into the future. Therefore, when mul- 
tiple descendants produce two branches reaching z — 0, we 
combine them into one galaxy if the split occurred within 
0.7 Gyr (half the dynamical time of the dark halo). If the 
two branches have remained separate longer than that, we 
assume that they will never rejoin, and instead separate the 
common progenitor into two galaxies, each one getting the 
descendant's fraction of the total mass of the SKID group. 
In these cases we always find that the SKID group has two 
parents with masses nearly identical to those of the two de- 
scendants; thus it appears that what has happened is that 
two galaxies have momentarily passed close enough to one 
another to be identified as one SKID group but do not have a 
strong physical interaction. The combining of unmerged de- 
scendants at z = means that we start off with 516 SKID 
groups that turn into 507 merger trees (galaxies). 

As is generally the case in simulations of this sort, ef- 
ficient cooling leads to galaxy masses at z = that are 
systematically larger than observational estimates. Figure 



2 We have checked that our trees remain essentially unchanged 
if we change this requirement to 8 or 32 particles. 




Figure 2. The figure shows the mass accretion histories for three example galaxies. In each row, the leftmost panel shows the masses 
of all of the galaxy's progenitors as a function of look backtime. When galaxies merge the lower mass progenitor is connected to the 
main progenitor by a vertical line. Vertical dotted lines mark rcdshifts z = 0.5, 1, and 2. The remaining panels show projected particle 
distributions at z = 0, 0.5, 1, and 2, in regions 0.63 comoving Mpc (top row) and 2.5 comoving Mpc (middle and bottom rows) on a side. 
Green points denote particles that end up in the galaxy at z = 0. Red and blue points mark other gas and star parciles, respectively. 



shows the z — cumulative mass function, <E>(M ga i), of 
this simulation in comparison t o the funct i onal fit to the 
z = mass function estimate of iBell et alJ l)2003bJ) . which 
uses galaxy luminosities and colors to infer stellar masses as- 
suming a "diet Salpeter" initial mass function (IMF). The 
cumulative mass function is just the number density of ob- 
jects with mass greater than M, or 

P oo 

$(M) = / (p(M')dM'. (1) 

J M 

The simulated galaxies are too massive at all number den- 
sities, with the discrepancy worsening at low and high mass 
values. The source of this discrepancy remains unclear even 
after many years of study . It is possible that it is p artly 
a numerical artifact (e.g. ISpringel fc Hern auisd l2 002t) . or 
the result of insuffi cient supernov a (e.g. lDekel fc Silkl ll986) 
or AGN (e.g. lOmma et, alJ l2004h feedback in the simula- 
tions, ortiie_^sjilt of inc orrectly assuming a universal stellar 
IMF. iKeres et ail (2005) argue that a physical mechanism 



that suppresses "hot" accretion from shock heated gas halos 
would substantially reduce galaxy masses, especially above 
M* . Such a mechanism might involve AGN feedback or con- 
duction, or it might sim ply emerge from a prope r treatment 
of multi- phase cooling IIMaller fc B ullock 2004), which can 
reduce the overall rate o f cooling onto massive galaxies (see 
iDekel fc BirnboimlEool for a recent review of all these pro- 
cesses) . 

Because we compute merger statistics in terms of galaxy 
mass ratios, we are hopeful that any solution to the global 
mass discrepancy will have only a modest effect on the 
galaxy merger statistics that we present here. There may, 
however, be some changes to merger rates because lower 
mass galaxies have longer dynamical friction timescales. 
We plan to investigate the importance of this effect in fu- 
ture work. The correction to galaxy masses needed to make 
the simulated mass function agree with the observations is 
shown in the bottom panel of Figure Q In the galaxy mass 
range that we consider, 2x 10 10 Af Q < M gal < 6x lO n M at 
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Figure 3. The number of mergers per galaxy as a function of 
redshift for the high mass sample. The three solid lines are, from 
top to bottom, -Rmcrg > 0.125, Rmcrg > 0.25, and i? mcrg > 0.5. 
The dashed lines show the maximum number of possible merg- 
ers including unresolved mergers. The diamonds mark the last bin 
where the Poisson error on the resolved mergers is larger than the 
possible contribution of unresolved mergers and the lines are con- 
tinued for one bin past this point. The dividing redshift between 
these bins and the next are z ros = 0.5, 1.0 and 2.6. 



z — 0, the correction is well approximated by a constant fac- 
tor of 2.75 (dotted line). Thus we expect that a physical or 
numerical fix to the mass discrepancy would have only mod- 
est impact on the merger rates at a given mass ratio. How- 
ever, the correction becomes more significant at low masses 
(a factor of 18 for the lowest mass bin), which suggests that 
we may significantly overestimate the number of mergers 
when one parent has a mass less than 1O 1O M0. We will use 
the galaxy masses from the simulations in our discussion, 
but when connecting to observations we note that to get 
the correct number density of galaxies these masses should 
be divided roughly by a factor of 2.75. An alternative way 
to make the connection to observed galaxy populations is to 
match the cumulative space density above mass thresholds 
in our simulation to the cumulative space de nsity computed 
from the observed galaxy luminosity function dBlanton et alJ 
2003), assuming that luminosity is a monotonic function 
of mass. The luminosities implied by this number density 
matching procedure are given in Tabled 



3 THE PROPERTIES OF GALAXY MERGER 
TREES 

Figure [5] illustrates the merger history of three simulated 
galaxies, selected to show a range of behaviors. In the left 
panels, lines show the baryonic masses (stars plus cold gas) 
of each galaxy's progenitors as a function of lookback time. 
In the second column, green points mark particles that are 
members of the galaxy at z = 0, and blue points mark stars 
in other galaxies. Red points mark the surrounding gas par- 
ticles, most of which are at high temperature. The remaining 
panels show the positions of these particles at z — 0.5, 1, and 
2, in regions of constant comoving size. Green points always 
represent particles that will be in the galaxy at z = 0. 



Figure 4. The number of mergers per galaxy as a function of 
redshift for the medium and low mass samples. The triangles 
denote the medium mass sample with the upper solid line be- 
ing -Rmcrg > 0.25 and the lower solid line being R me rg > 0.5, 
while the square marks the line for the low mass sample with 



Rn 



> 0.5. The dashed lines show the maximum number of 



possible mergers including unresolved mergers. The symbols mark 
the last bin where the Poisson error on the resolved mergers is 
larger than the possible contribution of unresolved mergers and 
the lines are continued for one bin past this point. The dividing 
redshift between these bins and the next are ^rcs = 0.5,0.6 and 
1.0 for the low mass sample and the medium mass sample with 



Rn 



> 0.25 and R me rg > 0.5, respectively. 



The topmost galaxy has three clearly identifiable pro- 
genitors at z = 0.5. It experiences two major mergers 
(roughly 1:3) in the final Gyr, and it still has a bimodal ap- 
pearance at z = 0. In fact, SKID splits the system into two 
components at z = 0, but they are a single component at an 
earlier output and are therefore counted as a single galaxy. 
The galaxy in the middle panels undergoes its last major 
merger at z ~ 1.2, and thereafter it grows by smooth accre- 
tion, doubling its mass by z = 0. The morphological type 
recipe that we adopt in [21 (similar to that in semi-analytic 
models) assigns the post-merger accretion to a disk compo- 
nent and therefore assigns this galaxy a bulge-to-total mass 
ratio of ~ 0.5. However, if s ome mechani s m suppresses hot 
gas accretion at late times llBin nev||20 04t iKatz et aljl2003l : 
iDekel fc BirnboimlEool iKeres et al J 120051) . then the bulge 
fraction would be higher (see 21 for further discussion). The 
galaxy in the bottom panel experiences many mergers, but 
all of them are minor; our morphological type recipe would 
therefore identify this as a disk-dominated system. 

In quantitative terms, we would like to know the func- 
tion where 



f(M gal ,z,R n 



OPrr. 



dtdM SBl dRn 



(2) 



the probability that a galaxy with a mass between M ga i and 
M ga i + dMgai at a redshift z was the result of the merg- 
ing of two galaxies with a mass ratio between i? mer g and 
fl merg + di? mcrg in the preceding time interval dt. We de- 
fine -Rmcrg < 1) i- e -> it i s t ne ratio of the baryonic mass of 
the smaller parent to that of the larger parent. Unfortu- 
nately, we do not resolve a large enough dynamic range to 
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Figure 5. The number of nearly equal mass mergers (i^merg ^ 
0.5) per galaxy, that occur for > z > 0.5 as a function of the 
galaxy mass M ga{ at z = 0. The number of mergers is strongly 
dependant on the galaxy mass. The dashed line shows a linear fit 
to log(Af„ a ;) with a slope of 0.0145. 

fully characterize the function ty, so instead we study sev- 
eral interesting "projections" of it to learn about its overall 
parameter dependence. 

With the morphological assignment recipe of ||H a 
galaxy's Hubble type is determined by the merger history 
of its "main" progenitor, the most massive progenitor at 
each output. If the Hubble types of the merging galaxies 
change the result of the merger, e.g. because bulges behave 
differently from disks, then the merger history of the minor 
progenitors also matters. However, this effect is likely to be 
of secondary importance, and in this paper we will not dis- 
tinguish between mergers of disk and bulge galaxies. TableQ 
lists the number of resolved main branch mergers (mergers 
involving a galaxy's main progenitor) and the total number 
of resolved mergers for all of our merger trees. Only for the 
high mass sample, M ga i > 6.4 x 1O 1O M0, do we resolve a sig- 
nificant number of off-main branch mergers. In most cases, 
therefore, we present results only for main branch mergers, 
but we occasionally compare the statistics of main branch 
mergers to those of all mergers for the high mass sample. 

3.1 Galaxy-averaged merger statistics 

Figures [3] and 2] present our main characterizations of 
$ (Mgai, z, -Rmcrg): the average number of mergers per galaxy 
per Gyr above Rmerg thresholds of 0.125, 0.25, and 0.5 for 
galaxies in the high mass sample (Fig. [3J and the medium 
and low mass samples (Fig.|lJ. The relation of diVmerg/di to 
*(M ga i,2,i? mcrg ) is 

= 1 ^Et s r7L tt * (Mi ' z ' jRm ^ )dEm ^ (3) 

where R m in is the mass ratio threshold, JV ga i is the number of 
simulated galaxies in the mass range M m i n to M max , dn/dM 
is the galaxy baryonic mass function, and n is the mean 
space density of galaxies in the mass range. We compute 
dN merg /dt by counting mergers in a redshift interval Az 



Figure 6. The fraction of mergers with parent mass ratios i? mC rg 
for galaxies in the high (diamonds), medium (triangles) and low 
(square) mass samples are shown for mergers in the redshift range 
< z < 0.5. The points are normalized by the width of the bin. 
The distribution of the high mass sample is well fit by a function 
oc i?morg plotted as a dotted line. For the medium and low mass 
samples we do not have enough data to constrain the shape of the 
distribution, so we do not know if this distribution is a function 
of galaxy mass or redshift. 

and dividing by the corresponding time interval At, with 
typical values At ftf lGyr. The lower mass ratio mergers 
(i.e., i?min = 1/8 or 1/4) can only be resolved for higher mass 
galaxies. Merger rates are substantially higher for massive 
galaxies — e.g., the average rate of -Rmerg > 0.5 mergers 
at z = 0.3 is 0.054 Gyr -1 for the high mass sample and 
0.018 Gyr -1 for both the medium and low mass samples. 

Our main limitation in computing these statistics is that 
we do not reliably resolve SKID groups with fewer than 64 
particles. Thus, when a galaxy's mass increases by 63 par- 
ticles we cannot tell without detailed examination whether 
this growth was the result of a merger or of smooth accre- 
tion. The directly calculated rates shown by solid lines in 
Figures [3] and [I] are therefore lower limits to the true rates. 
Dashed curves accompanying the solid curves show rates 
that include the maximum contribution of unresolved merg- 
ers, assigning all growth of up to 63 SPH particle masses that 
is not in resolved mergers to unresolved mergers. At low red- 
shift there are no unresolved mergers, but as we go back in 
redshift and the galaxy masses decrease, the number of pos- 
sible unresolved mergers increases. We stop each line when 
the number of possible unresolved mergers is greater than 
the Poisson uncertainty in the number of resolved mergers. 
Symbols mark the last output bin before this limiting red- 
shift, where we believe we resolve enough of the mergers to 
robustly study their properties. Since the assumption used 
for the dashed lines is extreme, we expect that the solid lines 
are generally closer to the true rates. 

For the high mass sample, the limiting redshifts for 
Rmerg > 0.125,0.25 and 0.5 are z rcs = 0.5,1.0 and 2.6, re- 
spectively. For medium masses, merger ratios of -Rmcrg > 
0.25 and 0.5 can be resolved to z las — 0.6 and 1.0, respec- 
tively. For the low mass sample, Rmerg > 0.5 can be resolved 
only to z rcs = 0.5. To compare all three samples, we must 
therefore restrict ourselves to Rmerg > 0.5 and z < 0.5. 
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Figure[!>]quantifies the mass dependence of merger rates 
seen in Figures[3]and2] To maximize statistics, we count all 
mergers with -Rmcrg > 0.5 occurring between z = and 
2 = 0.5. The mean number of mergers per galaxy in this 
interval follows a best-fit relation 



N„ 



= 0.0145 log 



10 n M, 



+ 0.0175, 



(5) 



where M ga i is the galaxy mass at z = 0. We can therefore 
anticipate a strong dependence of Hubble type frequencies 
on galaxy mass. We find an even steeper best-fit slope for 



Rn 



> 0.25, but with only two points it is hard to draw 



meaningful conclusions from this difference. 

Figure |S| shows the distribution of parent mass ratios 
for main branch mergers in the interval < z < 0.5. In the 
high mass sample, the number of mergers is proportional to 
^mcrg- There are not enough bins in the low and medium 
mass samples to reliably infer a dependence (or lack thereof) 
on galaxy mass. The —1.2 slope for the high mass sample 
implies that high mass ratio mergers dominate the merger 
growth rate, since the average rate at which a galaxy gains 
mass by mergers above a threshold i? m in, 



dM i 

— OC / *(i? me rg)i? n 



crdRn 



ocl-R" 



(6) 



has already reached half of its asymptotic value for R n 
0.42. 



3.2 Volume-averaged Merger Rates 

We now turn to the volume-averaged merger rate of galax- 
ies (the mean number of mergers per comoving Mpc 3 per 
Gyr), which can be estimated observationally by counting 
either recent merger remnants or close pairs that will merge 
in the near future, and dividing by an estimated merger 
timescale. The merger rate of all galaxies (not just the main 
progenitor branch) is shown in Figure [7J for mergers with 
Rmerg > 0.5. We plot the result for all branches because 
observations cannot determine which merger remnants will 
end up merging with other galaxies at future times. Dashed 
lines include the possible contribution of unresolved mergers. 
We continue curves past the limiting redshift z ICS adopted in 
Figures[!J]and[l]because the range of possible merger rates is 
still an interesting prediction, even though the contribution 
of unresolved mergers may be significant. Although we are 
likely to be missing some mergers above z rcs , we think that 
the resolved merger rates (solid lines) are probably closer 
to the true merger rate predictions than the dashed lines, 
which are upper limits based on extreme assumptions. 

The naked e rror bars i n Fig ure Q show the observa- 
tional estimates o flLin et alj i2004) derived from the DEEP2 
survey dDavis et alJ 120031) . for galaxies with luminosities 
— 19 < Mb < —21. These points lie below our model pre- 
dictions, especially at z > 1. The luminosity range roughly 
corresponds to the mass range of our high mass sample, ex- 
cept that these are the galaxy luminosities at the observed 
redshift, not at z = 0. If we restrict ourselves to galaxies 
that are in that mass range at the plotted redshift, then we 
get the dot-dash line with crosses in Figure Q which agrees 
fairly well with the observations at z < 1 but remains high 
by a factor ~ 2 at z > 1. Evolution of stellar mass-to-light 
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Figure 7. The number of mergers per Gyr per comoving Mpc is 
shown as a function of redshift. The high mass sample (diamonds) 
and medium mass sample (triangles) are shown for all resolved 
mergers (solid lines) and including the total possible number of 
unresolved mergers (dashed lines). We see that the contribution 
of unresolved mergers is negligible for z < 1.7 in the high mass 
sample and z < 1.0 in the medium mass sample. The merger rate 
when galaxies are selected by their mass at that redshift is shown 
as the dot-dash line with crosses. This is more relevant to what 
one sees obs ervationall y in a flux limited sample. For comparison 
the data of iLm" et alJ 120041 ) is shown as the naked error bars 
for a similar range in stellar mass. The error bars on the upper 
most points are representative of the errors for all the simulation 
points. 



ratios changes the correspondence between luminosity and 
stellar mass, and the range of mass-to-light ratios becomes 
larger at earlier epochs, so a full assessment of this mild 
discrepancy will require more detailed modeling of the stel- 
lar populations of the simulated galaxies, and more detailed 
replication of the observational procedures for estimating 
merger rates. The systematic uncertainties in matching the- 
oretical and observational samples are minimized if merger 
rates are given divided by the number density of objects 
being sampled, as in our Figures |3] and 2] When data are 
presented as in Figure [7] the number of mergers per galaxy 
is convolved with evolution of the number density of galax- 
ies and of their properties, making it more difficult to isolate 
the source of discrepancies. 

Figure |H] shows the distribution of the total number of 
mergers a galaxy undergoes during its history, for different 
mass samples and mass ratio thresholds. Figure legends list 
the average number of mergers per galaxy for each mass bin 
and parent mass ratio. Curves in each panel show a Pois- 
son distribution with the corresponding mean, and in all 
cases they describe the measured distribution with reason- 
able accuracy. This result suggests that each merger is an 
independent event; the average probability of mergers de- 
pends on galaxy mass, but in a given mass range there are 
not "merger heavy" or "merger light" galaxies beyond what 
is expected from Poisson statistics. 




number of mergers 

Figure 8. The number of mergers a galaxy undergoes for the 
low, medium and high mass samples and for parent mass ratios 
of 1 : 2 and 1 : 4. We see that all five are well fit by a Poisson 
distribution with a galaxy having a mean number of mergers of 
0.3, 0.2, 0.8, 0.7 and 1.3, respectively. 



3.3 Last Major Merger 

We now turn to another interesting characteristic of merger 
histories, the distribution of redshifts at which galaxies ex- 
perience their last major merger. Figure |§] plots the fraction 
of galaxies that have undergone a merger above a 1:2 or 1:4 
mass threshold since redshift z, in the three mass bins. In- 
cluding unresolved mergers has only a modest effect on this 
result and is most significant at higher redshifts (see dashed 
lines). We see that in all cases a sizeable fraction of galaxies 
have never undergone a merger since the time the galaxy 
first crossed our mass resolution threshold of 6.8 x 10 9 Mq. 
(Galaxies in the low mass sample have grown by a factor of 
3 — 5 since that time, while galaxies in the high mass sample 
have grown by at least a factor of nine.) In the medium mass 
sample, for example, only half of the galaxies experienced a 
merger with mass ratio larger than 1:4, and less than 30% 
had a 1:2 merger. 

The symbols in Figure [5] mark the median redshift of 
the most recent major merger for those galaxies that did 
experience a merger. About half of these mergers occur at 
z<l — 1.5. We can c ompare this resu lt to the predictions 
for dark matter halos. IWechslerl {2001), Fig. 137, shows the 
redshift distribution of the last 1:3 merger for dark matter 
halos in differing mass ranges. If we simply approximate the 
mass of the halo for each of the z — galaxies by multiplying 
by 2Q rn /£lb, we obtain median halo masses of 0.6, 1 and 2 x 
10 12 Mq for t he low, medium a nd high mass galaxy samples, 
respectively. 3 IWechslerl <200lft finds that 40% and 28% of 
halos with masses of 0.6 and 1.4 x 10 12 Mq, respectively, 
have not had a major merger (greater than 1:3) since z = 
4. For those halos that did have major mergers at z < 4, 
the median redshift of the most recent merger is z = 1 in 

3 For this halo mass range, we find in the simulation that roughly 
50% of the halo baryons end up in the central galaxy. 



Figure 9. The fraction of galaxies that have undergone a merger 
of parent mass ratio iJ mcrg > 0.5 (solid lines) or iJ mcr g > 0.25 
(dot-dash lines) as a function of redshift is shown for the three 
mass samples: high (diamonds), medium (triangles) and low 
(squares). We see that in all cases some galaxies never undergo a 
merger of that parent mass ratio. The corresponding dashed lines 
show the cumulative fraction including unresolved mergers. The 
median redshift of galaxies that do undergo a merger is marked 
by the symbols. 

both cases. Compared to the dark halos, more galaxies in 
our simulation are merger free, and the median r edshift of 
the last major merger is slightly higher. However. IWechsleil 
( 2001) uses a slightly different cosmological model (fi m = 
0.3, A = 0.7, erg = 1, and n = 1), which could plausibly 
account for differences of this magnitude. 

The simila rity of our galaxy merger rates to the 
IWechslerl ll200lll halo merger rates is broadly consistent with 
the idea that most major mergers of halos are quickly fol- 
lowed by mergers of their central galaxies. In future work we 
intend to perform a more detailed analysis of this point using 
halos identified in our simulation to investigate the relation- 
ship between dark halo and galaxy mergers. Unfortunately, a 
comparison of our results to semi-analytic models of galaxy 
formation is difficult, since most authors only present results 
that are meant to be compared directly with observations 
and do not present the intermediate results necessary for 
a comparison to the simulations or to other semi-analytic 
models. The only sem i-analytic predictio n with which we 
can compare is that of iBaugh et alJ (Il996l) . who show plots 
similar to Figure |5] In their model, nearly all galaxies have 
had a major merger since z ~ 2. However, they consider an 
Q m — 1 cosmology, which is undoubtedly the major cause of 
this large difference from our results. A comparison of Fig- 
urelSjto semi-analytic predictions using the same cosmology 
would be interesting. 



4 THE DISTRIBUTION OF HUBBLE TYPES 

We now turn to determining a galaxy's Hubble type 
from its merger history. Simulations have shown that 

collisions between disk galaxies typical l y resu lt in 

spheroidal systems llToomre fc Toomrelll972l : lBarne^ll988l : 
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Hernauistll99llBarnes fe Hernauistll996l : lNaab et alll999t 
Naab fc BurkertlbOOST) . We investigate a simple ansatz for 
the effect of mergers on morphology: if the mass ratio i? mC rg 
exceeds a threshold i? ma jor, we convert the entire baryonic 
mass of the merger remnant to a spheroid. 4 Mergers below 
the -Rmajor threshold have no effect. All mass accreted after 
the last major merger is assigned to the disk component. 

Figure ED shows the distribution of bulge-to-total mass 
ratios computed using this recipe and the simulated merger 
histories, for the three mass bins and -Rmajor thresholds of 
1:2 and 1:4. At low bulge-to-total ratio, unresolved merg- 
ers could have a significant impact. At the low end of our 
mass range we only resolve the later two thirds of a galaxies 
mass accretion history, so up to one-third of its mass could 
have had an unresolved merger before we start resolving it. 
We therefore show, in addition to our predictions using re- 
solved mergers, a prediction using the extreme assumption 
that each galaxy has the maximum possible contribution 
from unresolved mergers. In general, this maximal contri- 
bution simply shifts galaxies from the lowest bulge-to-total 
bin (< 0.1) to the next bin (0.1 — 0.2). Thus, our prediction 
for the fraction of truly bulgeless low mass galaxies can be 
significantly affected by unresolved mergers, but our other 
predictions cannot. 

For bulge-to-total ratios > 0.2, the predicted distri- 
butions in Figure 1101 are strikingly flat. In particular, the 
simulation produces relatively few systems that are close 
to 100% bulge, because after mergers galaxies continue to 
accrete gas and form new sta rs in a disk comp onent in 

the simulation. As discussed by Keres et al. (2005J; see also 
r— ii . 1 i— n i.. ,. _ 



Binnevll200l ICroton et all 12003: iDekel fc Birnboimll2004 . 

a mechanism that halts late-time accretion onto high m ass 
galaxies, such as A GN feedback j Binney fc Taborlll995l) or 
multi-phase cooling jMaller fc Bullockll200^ would improve 
the agreement between the simulations and observations by 
predicting lower masses and older stellar populations at the 
high end of the galaxy luminosity function. Reproducing the 
correct distribution of bulge-to-disk ratios could be a pow- 
erful test of such mechanisms, since the fraction of bulge- 
dominated systems will change depending on the conditions 
assumed for shutoff of gas accretion. 

As expected, the higher merger rates of massive galaxies 
lead to a greater frequency of high bulge-to-disk ratios in our 
highest mass bin. To summarize the mass dependence and 
compare to observations, we separate the galaxy populations 
into "early" and "late" subsets at a bulge mass fraction of 
0.5. For a major merger threshold of 1:2, the late-type frac- 
tion is about 85% in the low and intermediate mass bins but 
only 74% in the high mass bin. For a threshold of 1:4, the 
late-type fraction is 69% in the intermediate mass bin and 
46% in the high mass bin; we cannot resolve 1:4 mergers in 
the low mass bin. 

Figure [TTI compares these results to an estimate of the 

4 We have also investigated the recipe used in some semi-analytic 
models where for each merger of mass ratio -R mC rg > fimajor twice 
the mass of the less massive progenitor is converted into bulge, 
and this is done for all mergers along the galaxy's main branch. 
Final bulge-to-disk ratios for this "partial conversion" recipe are 
slightly lower than those for the "total conversion" recipe, but 
the qualitative results are very similar, so we do not present them 
separately. 
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Figure 10. The bulge-to-total fraction is shown for each mass 
sample and for two choices of what constitutes a major merger, 
Rmajor = 0-5 an d fimajor = 0.25. The solid line includes only 
resolved mergers and the dashed line includes all possible mergers. 
The shaded region shows the range that could exist given our 
finite resolution. Dividing early and late type galaxies at a bulge- 
to-total ratio of 0.5, we also indicate the late galaxy fraction, 
which is the same whether or not unresolved mergers are included. 



early- type fraction derived from Bell et al.'s (l2003bT) deter- 
mination of the early and late type stellar mass functions, us- 
ing th e 2MASS data set iSkrutskie et al.lll997tlMaller et alJ 
2005). Because of the offset between predicted and observed 
galaxy mass functions discussed in we divide the sim- 
ulated galaxy masses by 2.75 for this comparison; we are 
thus comparing populations of similar space density. For 
Rmajor =1:2, the predicted early-type fraction is substan- 
tially below the observational estimate. For 7? ma jor =1:4, 
the agreement is excellent, albeit with only two theoreti- 
cal data points. We have checked that ratios of 1:3 and 1:5 
are also in fairly good agreement. iBell et al. (2003b) deter- 
mine galaxy types from concentrations, and the principal 
systematic uncertainty in this comparison is that we do not 
know how well their definition corresponds to our theoret- 
ical criterion /bulge = 0.5. The full distribution of bulge- 
to-disk ratios would be a powerful test of our predictions 
and, more generally, a powerful diagnostic of the physical 
processes that produce galaxy bulges. Extracting unbiased 
estimates of this distribution for large, well defined galaxy 
samples is a c hallenging task ; but within reach of modern 
data sets ( e.g. iMcIntosh et alJEoolj) . 



5 CONCLUSIONS 

We have used an SPH simulation of a ACDM universe to 
calculate the statistics of galaxy merger histories. Our prin- 
cipal findings are as follows: 

1. The incidence of major mergers is much higher for 
high mass galaxies. At z = 0.3, we find a rate of 0.054 merg- 
ers per Gyr above a 1:2 mass threshold for galaxies with 
-Mgai > 6.4 x 1O 1O M0, but the rate for our lower mass galaxy 
samples is three times lower. 
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Figure 11. The fraction of galaxies that are early type as a 
function of the galaxies mass is shown for parent mass ratios of 
1:2 (solid line) and 1:4 (dashed line). T he crosses with erro r bars 
are the observed fraction derived from iBell et al . 1 2003b]). The 
masses of the simulated galaxies have been divided by a factor of 
2.75 to bring them into agreement with the observed stellar mass 
function. 



2. The distribution of merger mass ratios is N oc -Rmei-g 
in our high mass sample (the only one for which our statistics 
are sufficient to make the estimate). With this logarithmic 
slope, major mergers dominate the total mass growth from 
mergers, with half of the mass coming in mergers above a 
mass ratio thresho ld R mem = 0-4. Howeve r, as shown by 
iMurali et all (120021) and iKeres et alJ (120051) . even the high 
mass galaxies in our simulations gain most of their mass 
through smooth gas accretion rather than mergers with pre- 
existing galaxies. Our derived mass ratio distribution rein- 
forces their arguments that the dominance of smooth ac- 
cretion over mergers is not an artifact of limited numerical 
resolution. 

3. For each mass and mass ratio bin, the merger rate 
climbs rapidly with redshift, roughly doubling between z ~ 
and z ~ 0.6. The predicted merger rates are in reasonable 
agreement with observational estimates, but the systematic 
uncertainties in the comparison are large because of the dif- 
ficulty of matching simulated and observed galaxy popula- 
tions. 

4. In each bin of mass and mass ratio, the number of 
mergers per galaxy is Poisson distributed with respect to the 
mean number for the bin. This suggests that each merger is, 
effectively, an independent event. 

5. A substantial fraction of galaxies experience no major 
mergers after z = 1, and many have no major mergers over 
the entire time that they are resolved in our simulations. For 
example, even in our high mass sample, only 45% of galaxies 
have an i? mC r g >1:4 merger at z < 1, and only 20% have an 
i? merg >1:2 merger. For those galaxies that do experience 
major mergers, the median redshift of the most recent event 
is z ~ 1 — 1.5 in all mass and -Rmcrg bins. 

The high frequency of quiescent merger histories is per- 
haps our most important result. If major mergers are in- 
deed the primary mechanism of spheroid formation, then 
this large population of merger-free galaxies shows that the 
ACDM model can easily produce an acceptable fraction of 



late-type, nearly bulgeless systems. The predominance of 
quiescent merger histories is also good news for the galaxy 
angular momentum problem, since low redshift mergers are 
frequently responsible for removing angular momentum and 
shrinking disk sizes in simulations that focus on the forma- 
tion of individual galaxies. 

We have combined our merger histories with a sim- 
ple, semi-analytic style recipe for spheroid formation to es- 
timate bulge fractions of our simulated galaxies. With a 
bulge conversion threshold of 1:4, we find acceptable agree- 
ment with the observed trend of early-type galaxy fraction 
with galaxy mass, once we rescale the simulated galaxy 
masses to match the observed galaxy mass function. How- 
ever, since most galaxies have continuing gas accretion af- 
ter their last major merger, the fraction of truly bulge- 
dominated systems (i.e., /b u i ge >0.8) is small, even at high 
masses. This result suggests that some mechanism that shuts 
off gas accretion in high mass halos is needed to explain 
the observed distribution of galaxy morphologies. Many au- 
thors have already argued that a mechanism of this sort 
is required to reproduce the observed exponential cutoff 
of the galaxy luminosity functio n and the red c o lors o f 
the most ma s sive g a laxies (e . g., [ Kauffmann_j3t_alJ Jl993); 
Benson et alJ (120031) : iBinnevI (120041): [Croton et all 120051) : 
Dekel fc Birnboiml !2004l) : IKeres et al l J2005I) ). ~ 

Our simulation volume is too small to allow a detailed 
investigation of the morphology- density relation. However, 
in the one cluster mass halo (M v j r = 3 x 10 14 M©) halo that 
forms in our simulation, we do not find the high fraction 
of early-type galaxies seen in observed clusters. This dis- 
crepancy suggests that some mechanism other than major 
mergers must contribute to morphological transformations 
in the cluster enviro nment, such as ram pressure stripping 
JCunn fc GottHl972D or "harassment" by weak perturba- 
tions jMoore et alJll99r3) . 

We hope to extend this work in future studies that use 
simulations of larger dynamic range. These will allow us to 
investigate the correlation of merger histories with environ- 
ment and to better characterize the merger histories of lower 
mass galaxies. We will also investigate the relationship be- 
tween individual galaxy and dark halo merging histories, 
and we will check whether mechanisms that reduce galaxy 
baryon masses, such as galactic winds and suppression of 
hot accretion, also change the statistics of galaxy mergers. 
Comparison of our results to observations already suggests 
that the simplest picture of spheroid formation by mergers 
is not the full story. Future comparisons drawing on better 
simulations and more detailed observational analyses should 
yield much greater insight into the physics that determines 
galaxy morphologies. 
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